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Abstract 

The underdoped high- Tc materials are characterized by a competition between 
Cooper pairing and antiferromagnetic (AF) order. Important differences between 
the superconducting (SC) state of these materials and conventional superconductors 
include the d-wave pairing symmetry and a remarkable magnetic response to 
nonmagnetic perturbations, whereby droplets of spin-density wave (SDW) order can 
form around impurities and the cores of vortices. In a simple picture, whenever SC is 
suppressed locally, SDW order is nucleated. Within a mean-field theory of d-wave SC 
in an applied magnetic field including disorder and Hubbard correlations, we show in 
fact that the creation of SDW order is not simply due to suppression of the SC order 
parameter, but rather due to a correlation-induced splitting of the electronic bound 
state created by the perturbation. Since the bound state exists because of the sign 
change of the order parameter along quasiparticle trajectories, the induced SDW 
order is a direct consequence of the d-wave symmetry. Furthermore the formation 
of anti-phase domain walls is important for obtaining the correct temperature 
dependence of the induced magnetism as measured by neutron diffraction. 

* Corresponding author e-mail: markus.schmid@physik.uni-augsburg.de 
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A superconductor is characterized by a Bardeen-Cooper-Schrieffer (BCS) order param- 
eter Ak(R), where R is the center-of-mass coordinate of a Cooper pair of electrons with 
momenta (k, — k). The bulk ground state of such a system is homogeneous, but a spa- 
tial perturbation which breaks pairs, e.g. a magnetic impurity, may cause the suppression 
of Ak(R) locally. What is revealed when SC is suppressed is the electronic phase in the 
absence of Ak, a normal Fermi liquid. Thus the low-energy excitations near magnetic im- 
purities and in the vortex cores of conventional SC are essentially Landau quasiparticles 
trapped in bound states. The underdoped cuprates have been studied intensively in recent 
years in part because their proximity to the Mott insulator is thought to be responsible for 
many unusual properties, including possibly high-temperature SC itself. These systems are 
quite different from conventional SC, because when the pair amplitude is suppressed locally, 
e.g. by a vortex, a competing ordered state^ stabilized by the proximity to the Mott state 
appears to emerge instead of a normal metal. This state is characterized at low temperatures 
T by static local SDW order with an ordering wave vector Q near (7r,7r), an order which 
is not present in the state above Tp. This was reported first in elastic neutron scattering 
experiments^ on La2-xSr2:Cu04 (LSCO), with a correlation length of several hundred A, but 
has been confirmed in other underdoped cuprate materials as well."^"^ An enhancement of 
incommensurate static order was observed with increasing the applied magnetic field up to 
14T.^ Because the signal disappeared above T^, the magnetism was attributed to the vor- 
tices; indeed, scanning tunnelling microscopy (STM) measurements^ on Bi2Sr2CaCu208+5 
(BSSCO) were able to directly image unusual charge order near the vortex cores which is 
almost certainly related to the field-induced SDW detected by neutron scattering. 

Hints of magnetic ordering in the SC state had been detected earlier by ^SR experiments 
in zero field^'^ as a wedge-shaped extension of the "spin glass" phase into the SC dome of 
the temperature vs. doping phase diagram of cuprates (see Fig. la). Lake et al? reported 
that an incommensurate magnetic order similar to the field-induced state was observed in 
zero field, too. But although it also vanished at Tc, the ordered magnetic moment in zero 
field had a T dependence which was qualitatively different from the field-induced signal. The 
zero-field signal was attributed to disorder, but the relation between impurities and magnetic 
ordering remained unclear. Because strong magnetic fiuctuations with similar wavevector 
are reported at low but nonzero energies in inelastic neutron scattering experiments on these 
materials, e.g. on optimally doped LSCO samples exhibiting no spin-glass phase in zero field. 
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Figure 1: Phase diagrams and temperature dependence of magnetic order, a. 

Schematic temperature T vs. doping x phase diagram for cuprates. b. t/U dependence of the 
SC Tc and the disorder-induced magnetic transition temperature Tg for an impurity concentra- 
tion riimp = 10% with potential strength Vimp = 2.0t and a hole concentration x = 0.1. t is the 
nearest-neighbor hopping amplitude and U denotes the Hubbard on-site repulsion, c, d. Magnetic 
Fourier component |-/Vf(q)p at the ordering wavevector integrated around (vr,7r) vs. temperature 
for (c) different interaction strengths U (at fixed Vimp = 2.0i) and (d) different impurity potential 
strengths Vimp (at fixed U = 2.bt). 

it is frequently argued that impurities or vortices simply "pin" or "freeze" this fluctuating 
order. 

Describing such a phenomenon theoretically at the microscopic level is difficult due to 
the inhomogeneity of the interacting system, but it is important if one wishes to explore 
situations with strong disorder, where the correlations may no longer reflect the intrinsic 
spin dynamics of the pure system. Such an approach was proposed by the current authors 
in a model calculation for an inhomogeneous rf-wave SC with Hubbard-type correlations 
treated in mean field. In this model a single impurity creates, at sufficiently large Hubbard 
interaction U and impurity potential strength Vimp, a droplet of staggered magnetization 
with a size corresponding to the AF correlation length of the hypothetical pure system. ^^"^^ 
Such local impurity-induced magnetism has been studied extensively both theoretically and 
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experimentally, and was recently reviewed in Ref. 15. As was shown in Ref. 11, when 
these droplets come close enough to interact, there is a tendency to form incommensurate, 
phase-coherent Neel domains whose size is sufficient to explain the observations by Lake et 
al? in zero field. In addition, such a model explains the empirical observations that both 
increasing disorder^^ and undcrdoping^ enhance the SDW order. 

In this paper we investigate the origin of the "order by disorder" phenomenon described in 
Ref. 11, as well as the T evolution of the disordered magnetic state in apphed magnetic field. 
In the case of the field-induced SDW, an apparently very natural approach to the problem 
was developed by Demler et al.,^"^ who constructed a Ginzburg-Landau (GL) theory for 
competing SDW and SC order in a magnetic field. This order-parameter phenomenology 
describes correctly the reduction of condensation energy in the vortex phase of the pure SC, 
and leads to a phase diagram qualitatively consistent with experiments.^ The GL approach, 
however, ignores the energies of the quasiparticles moving in the inhomogeneous state which 
also can crucially affect the competition between SDW and SC order in these materials at 
low T, as we show here. 

In a d-wave SC without AF correlations, a bound state of an isolated vortex is found 
at zero energy^^ due to the sign change of the order parameter on quasiparticlc trajectories 
through the vortex core. On the other hand, solutions of the Bogoliubov-de Gennes (BdG) 
equations describing coexisting d-wave SC and SDW order^^"^^ show that this resonance is 
split by the formation of the SDW; that is, the system can lower the energy of the nearly 
bound quasiparticles by moving them below the Fermi energy. This finding is consistent with 
the STM experiments in the Abrikosov state of YBa2Cu307_5 (YBCO) by Maggio-Aprile 
et al!^^ and of BSCCO by Pan et al.,^'^ which observed split peaks in the vortex cores. A 
similar bound state is associated with non-magnetic impurities such as Zn in BSCCO and 
has been also imaged by STM.^^ It is therefore important to explore the role of quasiparticle 
bound states and their coupling to the SDW order to identify the origin of both types of 
induced local AF in the SC state. A more detailed understanding of field-induced order 
is also highly relevant for the interpretation of the quantum oscillations observed in recent 
transport experiment in high magnetic fields. ^^'^"^ These oscillations are possibly due to the 
formation of Fermi surface pockets as a consequence of SDW ordering and the concomitant 
reconstructed bandstructure. In addition, finding ways to understand the effects of disorder 
is crucial in order to reveal the intrinsic AF correlations present in the underdoped part of 



4 



the cuprate phase diagram, not least because the magnetic fluctuations at higher energies 
may be responsible for SC itself. 

The inhomogeneous mean-fleld theory presented here for electrons hopping on a square 
lattice with a rf-wave pairing potential and subject to a Hubbard on-site repulsion U, repro- 
duces the essential aspects of the field-induced spin-glass phase shown schematically in Fig. 
la. The primary purpose of this analysis is to model the inhomogeneous SC state and does 
not intend to describe the Mott transition to an insulating state at half-fiUing; therefore 
the doping dependence may not be directly compared to experiment. On the other hand, 
if doping is assumed to be correlated with the ratio of bandwidth to local Coulomb repul- 
sion t/U in the model, a phase diagram very much like the one found in various cuprate 
materials is obtained, as shown in Figs. Ib-d. We consider this a reasonable qualitative 
approach, since changes in the Fermi surface of these materials reported by angle resolved 
photoemission spectroscopy (ARPES) over the "spin glass" doping range are small, and 
it is therefore plausible that the primary effect on the electronic structure is due to the cor- 
relation induced band narrowing, as discussed in Ref. 11. The phase diagram we obtain in 
Fig. lb is thus comparable to the T — x phase diagram shown in Fig. la. The magnetically 
ordered phase can be enhanced both by the increase of the correlation strength or stronger 
disorder potentials, as shown in Figs. lc,d. 

The problem studied here involves several length scales, in particular the inter-impurity 
separation, the inter-vortex separation, the SC coherence and the magnetic correlations 
lengths, which can be difficult to disentangle. We obtain results which reproduce well the 
quahtative aspects of the experiment by Lake et al.,^ but show that some features depend 
on nonuniversal aspects of disorder, in particular the process of domain wall nucleation. 
While disorder- and magnetic-field induced SDW order both add to the ordered moment, 
the interference of disorder and magnetic- field effects is quantitatively significant. The do- 
main wall formation proves responsible for the distinct T dependences of the field- and the 
disorder-induced magnetization. An intriguing aspect of the present theory is that it also 
includes a crossover from magnetic droplets to filamentary stripe-like structures in selected 
regimes of hole densities and impurity concentrations. The model therefore offers a concrete 
route to describe the physics of the pinning of stripe correlations in the SC state. This 
insight may prove relevant for many experiments in the underdoped cuprates which have 
been attributed to stripes. 
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The basis for our model analysis is the BCS pairing Hamiltonian for a d-wscve SC with 
orbital coupling to an applied magnetic field, to which we add site-centered chemical dis- 
order and a local Hubbard repulsion; the latter is treated in an unrestricted Hartree-Fock 
approximation: 

ija ia (ij) 

+ ^ E ((^^)^^ - ("^^^ + E ^r'4c... (1) 

i ia 

Here, creates an electron on a square lattice site i with spin a =t, i- The hopping matrix 
elements between nearest and next-nearest neighbor sites i and j are denoted by tij = t 
and tij — t', respectively. An electron moving in the magnetic field from site j to i acquires 
additionally the Peierls phase cpij — (tt/^o) /r ' A(r) • dr, where $o = hc/{2e) denotes the 
superconducting flux quantum and A(r) = B{0,x) is the vector potential of the magnetic 
field in the Landau gauge. The chemical potential fi is adjusted to fix the average charge 
density n = Si(^j) = 1— ^c, where x is the hole concentration and N is the total number of 
lattice sites; in the following we will focus on a; = 0.1. The magnitude of the ci-wave pairing 
amplitude is determined by the strength of an attractive nearest-neighbor interaction 
Vd- The non-magnetic impurity potential V^^^ is described by a set of pointlike scatterers 
at random positions, and all fields, i.e. the pairing amplitude A^-, the local charge density 
(ni), and the local magnetization (af) are calculated self-consistently from the solutions of 
the associated BdG equations. (For further details of the numerical method we refer to the 
Supplementary Information.) 

We start with a single non-magnetic impurity in a d-wave SC at t/ = 0. The fingerprint 
of the induced virtual bound state is a single near zero-energy peak in the local density of 
states (LDOS) (see Fig. 2a). In this situation there is no magnetization induced by the 
impurity. Increasing the on-site Coulomb repulsion beyond a critical value Uc-, a staggered 
magnetization emerges in the neighborhood of the impurity (Fig. 2b). The two-sublattice 
nature of the magnetic pattern, in conjunction with the spatial extent of the impurity- 
induced resonance leads to a spin-dependent splitting of the near-zero-energy peak in the 
LDOS; one resonant state with a selected spin direction is thereby shifted below and the 
other above the Fermi energy, as is most clearly seen at the nearest-neighbor sites of the 
impurity where the resonant state has the largest weight (see Fig. 2a). The spin-dependent 
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Figure 2: Impurity- and field-induced magnetization, a, c. Local density of states (LDOS) 
at a nearest-neighbor site of a single impurity (a) and in the center of one vortex (c). In both cases, 
above a critical Uc, local SDW order is induced, concomitantly the zero-energy peak in the absence 
of the Hubbard repulsion (blue curve) splits spin-dependently (red curve). The dashed curves show 
the clean LDOS far from the perturbation, b, d. Real-space patterns of the magnetization (af) 
in units of hb on a 38 x 38 lattice for U = 2.2t > Uc at low temperature T = 0.025t. b shows 
the magnetization nucleated by a strong impurity (Vimp = 60t) located at the center. In d two 
superconducting flux quanta ^ = 2$o thread an impurity- free d-wave SC. 



splitting reduces the bound-state energy of a spin-up or -down state, and therefore stabihzes 
a local droplet of staggered magnetization which carries a total spin 1/2- moment surrounding 
the impurity. The splitting of the resonance peak can therefore be viewed as the origin of 
the impurity-induced magnetization. It is important to note that the splitting of the bound 
state is not due to the suppression of the rf-wave order parameter near the impurity. As 
we have verified, the order parameter can be artificially held constant in the solution of the 
BdG equations, and a nearly identical result is obtained. 

In the presence of a thermodynamically finite density of impurities, we recover at T = 
the results of Ref. 11, i.e. the creation of a defected but magnetically ordered state 
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as defined by strong peaks in the Fourier transform of the local magnetization M(q) at 
incommensurate wavevectors q near (tt, tt). The effect of temperature is now naturally 
included in the theory via T-dependent occupation probabilities of Bogoliubov quasiparticle 
states, as discussed in detail in the Supplementary material. In Fig. lb we show the extent 
of the quasi-ordered phase, labelled "spin glass" , which expands as correlations increase, and 
disappears at a critical value olt/U which depends on the strengths of the pairing interaction 
and the impurity potential. The intensity of the incommensurate magnetic Bragg peaks is 
shown in Figs. Ic and Id as functions of T for fixed t/U and fixed impurity potential 
strength, respectively. From Fig. Ic it becomes evident that the magnetic ordering, or 
"glass transition" temperature Tg can be smaller, equal to, or larger than T^. depending on 
U . In Fig. Id we show how increasing the impurity potential Vj^p can increase both the 
amplitude of the disorder-induced SDW and Tg itself. These results are consistent with the 
empirical observation that the size of the spin-glass phase is not universal, and in particular 
the critical doping beyond which magnetic order is no longer observed varies considerably 
between intrinsically disordered cuprates like LSCO and BSCO, and clean materials like 
YBCO. 

The magnetization induced by an orbital magnetic field can be traced to the same micro- 
scopic origin as the impurity-induced magnetization. ^° Above a critical Uc a staggered spin 
pattern is nucleated in the vortex cores with a spatial extent reaching beyond the size of 
a vortex core (see Fig. 2d), as observed in experiment.^ For the parameter set chosen, the 
core radius estimated from the area where the order parameter is suppressed is about one 
lattice spacing. The LDOS in the vortex center reveals that the origin of the field-induced 
magnetization is tied to the spin-dependent splitting of the Andreev bound state in the 
vortex core (Fig. 2c). The conjecture that the field-induced magnetization indeed appears 
simultaneously with the peak splitting in the LDOS is explicitly verified in Fig. 3. For the 
unmagnetized vortex at T = 0.1 75t a single Andreev bound-state peak exists at zero energy. 
With decreasing T the vortices nucleate a staggered spin pattern precisely at the T where 
the zero-energy peak in the LDOS sphts. With further coohng the peak splitting grows, 
more spectral weight is shifted below the Fermi energy, and the magnetization is enhanced. 

The natural next step is to consider a finite density of non-magnetic impurities in the 
presence of an external magnetic field and to compare it to zero-field results. Specifically for 
the modeUing of LSCO, we assume in the following that the Sr ions are the primary source 
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Figure 3: Temperature dependent peak splitting in a magnetic field. LDOS at the vortex 
center for three different temperatures below in the absence of impurities. Below the critical 
temperature Tg = 0.175i, the vortices magnetize and simultaneously the Andreev bound-state peak 
splits. The black dashed curve shows the LDOS far away from the vortex. 

of disorder, such that n^mp = a;, where n^mp denotes the impurity concentration. These 
systems are in the strongly disordered regime where the AF correlation length (droplet size) 
is comparable to the average distance between the dopants, such that the disorder is far 
from the one-impurity limit. Since the Sr dopants are removed from (but close to) the Cu02 
planes, we model them as weak scatterers with Vim-p = l-3t. Figs. 4a,b show the magnetic 
structure factor S'(q) at a fixed temperature far below Tc in zero and finite magnetic field 
averaged over ten different impurity configurations. As in Figs. lc,d, ^(q) is approximated 
by |M(q)p with M(q) = ;^ Xlj (^^^ Supplementary information for details). 

The magnetic signal in the structure factor appears at the incommensurate wavevectors 
q = (tt, TT ± 5) and q = (vr ± 5, tt) (see Fig. 4a). The magnitude of the incommensurability 
6 however varies for distinct impurity configurations randomly selected for 22 x 22 lattice 
systems. For the weak magnetic signal in zero field the averaging over different impurity 
configurations is therefore imperative but computationally demanding. Applying an external 
magnetic field strongly enhances the magnetization and reinforces incommensurate peaks 
at unambiguously selected wavevectors which are robust against variations in the impurity 
configurations (See Fig. 4b). 

Remarkably, the temperature dependence of the structure factor (see Fig. 4c) closely 
resembles the neutron scattering data on LSCO by Lake et al? For the results shown in Fig. 
4 we have chosen a parameter set where the staggered magnetization in zero field has its 
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Figure 4: Averaged magnetic structure factor for a many-impurity system. a,b. 

Intensity plot of the magnetic structure factor around (vr, vr) at T = 0.025t in zero magnetic field 
(a) and at finite field (b). The structure factor data were averaged over ten different impurity 
configurations. For the used system size of 22 x 22 lattice sites a magnetic flux of 2<&o corresponds 
to a strong magnetic field with H = 59T. The impurity concentration Uimp = x is fixed to 10% 
{Vimp = l-3t) and U = 2.9t > Uc- c. T-dependence of the peak intensity integrated around 
(vr, it) in zero-field and at finite field with the finite density of impurities riimp = x (blue and red 
curve, respecctively); for the data with $ = 2<I>o and riimp = x the zero-field data were subtracted. 
For comparison also the structure factor in a clean system is included for the same magnetic- 
field strength (purple curve). |M(q)p (integrated) translates directly to the ordered spin moment 
squared in units of fiB per Cu^^. 

onset at a temperature Tg indistinguishable from T^. This reflects a situation where upon 
cooling through Tc the localized bound states inside the d-wave energy gap emerge and 
immediately split in the self-stabilizing staggered magnetic pattern. Towards lower T the 
magnetic structure factor rises in a markedly different way in zero and in finite field. While 
the field-induced part of the magnetic signal has a negative curvature, the zero- field magnetic 
structure factor increases approximately linearly upon cooling. The two mechanisms of 
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Figure 5: Anti-phase domain walls, a, b. Real-space image of the staggered magnetization 
at T = 0.025t induced by three non-magnetic impurities. In impurity configuration a anti-phase 
domain walls appear vertically. In configuration b the staggered magnetization induced by the 
three impurities adjusts to a uniform AF domain around them. c. Temperature dependence of 
the integrated magnetic structure factor for the impurity configurations a (blue curve) and b (red 
curve), respectively. 



impurity- and field-induced SDW do not simply cooperate additively; the field-induced part 
of the magnetization is twice as large in the presense of impurities as compared to the field- 
induced SDW in the clean system (see Fig. 4c). The zero-field increase of the magnetization 
in the inhomogeneous SC state originates from the merging of AF patches nucleated by the 
individual impurities. Without impurities the increasing field-induced magnetization with 
cooling results from the growth of the magnetized regions around the well-separated vortex 
cores. Thus, both our zero- and finite-field results for the finite density of non-magnetic 
impurities riimp = x closely follow form of the T dependence of the the neutron scattering 
data for underdoped LSCO.^ Still, due to computational restrictions we are not yet able to 
access the low magnetic-field strengths to allow for a direct comparison with experiment. 



11 



An important remaining question is why the T dependences of the magnetization are 
different in zero and in finite field. A hint is provided by the observation that the T depen- 
dence of the magnetic structure factor for the impurity-free field-induced magnetization and 
also for just two single impurities in zero field has a negative curvature. In both cases the 
induced staggered magnetization patterns around each impurity or each vortex, respectively, 
adjust their individual two-sublattice spin structures in phase and thereby avoid any domain 
walls. For three nearby impurities, however, it proves already difficult to find a specific 
configuration where anti-phase domain walls are absent. In Figs. 5a,b we compare the stag- 
gered magnetization of two 3-impurity configurations with distinctly different domain-wall 
patterns. Remarkably, placing the three impurities on the same sublattice to form a right- 
angled equilateral triangle as in Fig. 5a generates a sequence of vertical anti-phase domain 
walls. If instead the impurities are configured in an acute equilateral triangle as in Fig. 5b 
a simply connected AF island forms around them. As Fig. 5c shows, with decreasing T the 
magnetic signal evolves differently for each impurity configuration. Intriguingly, |M(q)p 
rises almost linearly for the configuration with vertical anti-phase domain walls while it has 
a negative curvature for the single domain island. These examples, and others we have 
investigated, suggest that the linear low-T rise of the magnetic signal for a finite density of 
impurities in zero field originates from the anti-phase domain walls which are always present 
in the randomly generated impurity configurations. For the field-induced magnetization it 
is the larger distance between the magnetized vortices which prevents the occurrence of 
domain walls and therefore alters the T dependence of |M(q)p. 

All the results presented above focused on static disorder- and field-induced SDW, but 
inelastic neutron scattering experiments have shown that in the SC state the spin excitations 
at finite energy have almost the same distribution of spectral weight in q as the frozen 
magnetic order. For very low doping in the normal spin glass phase above Tc, the neutron 
intensity pattern is rotated by 45° and the connection to the spin correlations discussed here 
is less obvious. In fact, the utility of our model of choice is questionable for the description 
of the normal state where Fermi liquid concepts may not even be applicable. Nevertheless 
in the SC state we have provided a concrete foundation for the freezing of fluctuating spin 
correlations by disorder and magnetic flcld on the same footing; in particular the role of the 
quasiparticle bound states in the formation of the magnetic order has been highlighted. 

The new picture that emerges is complementary to the global competition between SC 
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h. U = 4:M,V^mp = IM, doping x = 15%, n'"^P = 7.5%, and Vd = 2.0t. In (a) and (b) the 
temperature is T = 0.025t. 



and AF phases in the sense that SC and disorder may significantly enhance SDW order in the 
underdoped regime. The ci-wave pairing of the SC condensate is crucial for this generation 
of local magnetism, as we have shown. Support for this cooperative effect between SDW 
and SC comes not only from the onset of the elastic magnetic neutron signal at but also 
from Zn-substituted optimally doped LSCO. There it is found by /iSR that 2% Zn induces 
a magnetic signal, but 3% Zn is found to eliminate it, but also destroys superconductivity^^; 
within the context of the current theory, this effect is understood not as a consequence of 
spin dilution'^^, but rather due to the destruction of the SC phase and thereby its ability to 
generate (or enhance) magnetic order. 

Finally we show that a qualitatively different kind of inhomogeneous textures may also be 
stabilized within the present weak-coupling approach. Figure 6a shows the typical charge- 
density profile in zero field for a parameter set used above to explore the onset of static AF. 
As expected, at and near the impurity sites the electron density is reduced, and in these areas 
the local SDW patches nucleate. With increasing repulsion U and for larger hole densities 
exceeding the impurity cocncentration, the inhomogeneous spin and charge patterns change 
qualitatively, and the impurity-centered patches with reduced electron density evolve into 
hole-rich filamentary structures (see Fig. 6b). In this still SC state the filaments constitute 
snake-like paths through an SDW background with an average density of almost one electron 
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per site. These textures provide a link to the study of disordered (quenched) stripes similar 
to those discussed recently within various GL models. ^^"^^ Therefore, depending upon the 
correlation strength and the details of the disorder, the magnetic ordering temperature Tg 
can vary significantly, and the ordering itself can be droplet-like or filamentary-like. This may 
explain much of the variability of neutron and IJ,SR experiments on different cuprates. Many 
interesting open questions remain to be addressed in future work, including the possibility 
of nematic instabilities in the presence of weak symmetry breaking fields and the transfer of 
spin-fluctuation spectral weight to flnite energies in samples where magnetism is not frozen. 

APPENDIX A: NUMERCIAL METHOD 

In order to investigate disorder- and field-induced magnetic order in rf-wave supercon- 
ductors we self-consistently solve the Bogoliubov - de Gennes (BdG) equations on a square 
lattice for the Hamiltonian 

ija ia (ij) 

+ I E - + E ^"^4aCi., (Al) 

i in 

where the hopping amplitude between nearest neighbor and next-nearest neighbor sites i 
and 2 is described by tij — t and tij — t' — — 0.4i, respectively. An orbtial magnetic field is 
represented by the Peierls phase factor ip^j — (tt/^q) /r ' A(r) • dr, while $o — hc/{2e) is the 
superconducting flux quantum and A(r) = B{0,x) is the vector potential of the magnetic 
field in the Landau gauge. The d-wave pairing potential is defined on two nearest neighbor 
sites i and j by 

\j = -Vd{cjiCi^) = Aji, (A2) 

where Vd is the attractive pairing interaction strength, which we set to = l.SAt throughout 
the paper. We then define a gauge invariant d-wave order parameter on each lattice site i 

= ^ (A?.+, + At_. - Aj,^. - AJ,_^) , (A3) 

where Af j — exp[— i(7r/$o) J^^ A(r) • dr\. The chemical potential is adjusted to fix the 
average charge density n—j^ while the electron number operator for spin cr at site i 
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is given by riia — cj^c^^, and the local charge-density operator by rii — respectively. 
Sf — \(tI — \ (nj-f — riii) is the z-component of the spin-operator at site i. 
The Bogoliubov transformation 

= X] {^inalnc + Kncll-^ > (A4) 
n 

diagonalizes the Hamiltonian in equation (Al), which thereby takes the form 



is the ground-state energy and 7^^^ creats an elementary fermionic Bogoliubov quasipar- 
ticle excitation with quantum number n, spin o", and energy Ena > 0. Calculation of the 
commutators of H from equation (A5) with the electron operators Qo- leads to a Schrodinger- 
like set of BdG equations 



E 



and 



E 



A* -H~ 





(A6) 




P I *"T 



(A7) 



with 



(A8) 



As we only search for solutions with positive E^a-, it is sufficient to solve the following single 
matrix equation 





(A9) 



This is because the solutions for En > are obviously identical to the solutions for the 
(positive) eigenvalues Eni equation (A6) 
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Figure 7: Division of the lattice into identical supercells. Lattice sites belonging to the same 
supercell are connected via solid lines, while dashed lines link sites of different supercells. Red 
lattice sites simulate a possibly existing disorder. 



while for En < the following relation holds 

and Eni = -E^ > 0. (All) 



En<0: 1 














[Vin J 



Since the solutions of equations (A6) and (A7) can be mapped on to those of the BdG 
equation (A9), we diagonalize equation (A9) to obtain the pairing potential Ajj, the charge 
density (nj), and the local magnetization (a^) self-consistently from 

^ij = -J 5Z i^i^^jn + UjnV*n) tauh (^^j , (A12) 
(n,t) = 5^|m,„|V(K), (A13) 

n 

K) = $^|t;.„P(l-/(i?„)), (AM) 

n 

(^n = ('^.t) - (^^i), (A15) 

where f{En) = (1 -\- e^"^''^^)~^ is the Fermi distribution function and T is the temperature. 
Sums over n run over positive and negative energies En- 

To maximize the size of the system for which equation (A9) can be diagonalized numer- 
ically, we take advantage of the magnetic translation symmetry of our model Hamiltonian 
(Al) by dividing the lattice into x My identical supercells each with x A^^^ sites (see 
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Supp. Fig. 7). 20,35,36 define the following magnetic translation operator^''' 

Tr = exp (-i R • (k + ^ A)) , (A16) 

where R is the translation vector and 7r translates any lattice vector r to the position 
r + R. Because [H, T^] — 0, it is possible to block diagonahze the Hamiltonian H in 
equation (Al) using the eigenstates of T^. This reduces the eigenvalue problem (A9) of 
dimension 2MxNx x 2MyNy to x My eigenvalue equations of dimension 2Nx x 2Ny. 
Applying the magnetic Bloch theorem 

block diagonalizes the BdG equations (A9), where k = 27r{j^^, j^^), Unk{ri) — Uink and 
■ynk(i"i) = Vink- Thus wc havc to solve the following 2Nx x 2Ny matrix equation for each k 
value 



E 

j 



^iJ+(k) A,,(k) 
^A^(k) -Hr.*{]^) 





(Ujnk^ 


— Enk 


Uink 1 


) 










\Vink j 



(A18) 



where 



^ nk 

K) = E M'fiEnk), (A20) 



^^^i) = jTw E i^-xr(i - /(Kk)). (A21) 



^ ^ nk 



i^ij and Ay are only k dependent, if i and j belong to different supercells. Then the back- 
mapping (see Supp. Fig. 8) leads to an additional phase for the u's and v's according 
to (A17), which is assigned to the matrix elements tjj(k) and Ajj(k). To make sure that 
two magnetic translations commute, we have to choose the magnetic field such that its flux 
through every supercell is a multiple of 2$o-^°'^^ Hence 2$o provides a lower boundary for 
the magnetic flux threading each supercell, which corresponds for a supercell enclosing an 
area of e.g. 22a x 22a to a magnetic field of about 59 T (we assumed a typical value for the 
in-plane lattice constant a in the cuprates of about a — 3.8A). 
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Figure 8: Hopping between supercells. A particle, which hops from the left supercell into the 
right supercell, is mapped back to the left supercell through the translation vector R. As a result 
the wave functions u and v obtain an additional phase given by the magnetic Bloch theorem (A17). 

In order to make contact with neutron scattering experiments, we evaluate the magnetic 
structure factor S'(q). In homogeneous systems it is defined as 

^(q) = ^E(^^'^o)^"''''^''"'"°^- (A22) 

i 

We approximate the spin-spin correlation function by the following factorization 

(^X)-(^nK)- (A23) 

Because the system which we are interested in is in general inhomogeneous, we have to sum 
over all lattice sites. Hence we find the expression 

= (A24) 

This approximation of the magnetic structure factor is identical to the Fourier transform of 
the magnetization squared. 

In Supp. Fig. 9 results for (nj), Af, and (a^) are shown in zero field (left column) 
and in finite magnetic field (right column) for a typical impurity configuration. One can 
identify the location of the impurities by the point-like suppression of the charge density 
(top row). While the d-wave order parameter is nearly homogeneous in the zero-field case 
(see Supp. Fig. 9b), one can clearly spot the positions of the two vortices where Af is 
suppressed to zero in Supp. Fig. 9e. In comparison to the zero-field case, a finite orbital 
magnetic field leads to an additional reduction of the order parameter over the entire lattice. 

18 



$ = $ = 2$ 




5 10 15 20 5 10 15 20 



Figure 9: Switching on an orbital magnetic field, a-c. Zero-field data. d-f. Finite-field 
data (<I> = 2$o)- d show the charge density (rij), b, e the d-wave order parameter A^, and c, 
f the magnetization (af) in real space. The same set of parameters is used here as in the rest of 
the paper, i.e. x = 10% = nimp,U = 2.9t,Vimp = l.3t. These data were obtained at the lowest 
temperature T = 0.025t we considered throughout the paper. Note the different scales in c and f. 



19 



Finally, for the parameters used here, the zero field case already contains impurity-induced 
antiferromagnetic order (see Supp. Fig. 9c), which is significantly enhanced by switching 
on a magnetic field. The magnetization peaks near the vortex cores, but due to the fact 
that strong type-II superconductors are penetrated by the field much beyond the cores, 
the magnetization is also enhanced in regions far away from the vortices, where the order 
parameter is nearly homogeneous. The SDW emerges due to the splitting of the Andreev 
bound state as explained in greater detail in the main body of the paper. 
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